An Exact Solution to the Quantized Electromagnetic Field in the de Sitter Space 
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In this work we investigate the quantum theory of light propagating in a de Sitter space-time. 
To do so, we use the method of dynamic invariants to obtain the solution of the time-dependent 
Schrodinger equation. The quantum behavior of the electromagnetic field in this background is 
analyzed. The quantum dispersion, the uncertainty product and the correlation between coordinate 
and momentum is evaluated. 
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INTRODUCTION 

The study of fields in a curved background has been 
considered along the last decades [l| . The quantization of 
gravity can be carried out by using the String Theory |2|- 
4]. However, in cosmological scales, gravity can be con- 
sidered as a classical theory and the fields as propagating 
waves in the background. Time-dependent backgrounds 
are used to describe many physical systems yielding in- 
teresting results. For instance, in the study of black hole 
evaporation [5| and the Unruh and Casimir effects [y, [D] . 
They are also very useful to describe the dynamical evo- 
lution of the universe, where the production of particles 
in cosmological spacetimes has been investigated [8rtL0( . 

The quantum effects of a massive scalar field in the 
de Sitter spacetime from a Schrodinger-picture point of 
view has been invstigated in Ref. [ll| . The authors 
used an exact linear invariant and the Lewis and Riesen- 
feld method [12l j to obtain the corresponding Schrodinger 
states in terms of solutions of a second order ordinary 
differential equation. By constructing Gaussian wave 
packet states, they calculated the quantum dispersions, 
quantum correlations, and the uncertainty product for 
each mode of the quantized scalar field. Scalar fields 
have also been investigated in Ref. [13f, where the au- 
thors construct the coherent states and established the 
existence of squeezed states. 

The quantization of the electromagnetic field in the 
de Sitter spacetime could provide a better understand- 
ing of (time-dependent) cosmological models, e.g., the 
present distribution of the fluctuation spectrum of the 
electromagnetic fields can be traced back to the primor- 
dial spectrum of cosmological fluctuations. 

In the context of the invariant method Pedrosa et 
al. [lj] investigated the light propagation through time- 
dependent dielectric linear media in the absence of free 
charges and in a curved spacetime from a classical and a 
quantum point of view. They found that the light prop- 
agation present a remarkable similarity in both cases. 
From the classical study, they showed that the ampli- 
tude q of the potential vector ( A ) is written in terms of 
the first- and second-order Bessel function for a propaga- 



tion through a dielectric media with e(t) — eoe at (a > 0). 
From the quantum perspective, they quantized the elec- 
tromagnetic waves propagating in a material medium 
with e(t) — £oe Qt , and q is written in terms of har- 
monic functions of exponentially decaying argument for 
a propagation in the de Sitter spacetime with the met- 
ric a(t) = e Hot , where Hq is the Hubble constant. They 
wrote the solutions of the Schrodinger equation in terms 
of p, which is solution of the Milne-Pinney equation 
[la, [iff . As a result they claim that for a problem with 
curved spacetime the electromagnetic field quantization 
can also be carried out following the same steps used for 
the time- varying dielectric case. Moreover, as in the clas- 
sical case, the quantum behavior of the electromagnetic 
field is purely oscillatory in the conformal time in the de 
Sitter spacetime. 

The study of de Sitter spacetime becomes important 
if one considers a ACDM model for the cosmos. In the 
present era, the universe can approximately be described 
by a de Sitter spacetime, where the matter decay with the 
volume and the cosmological constant is kept constant. 
For t ^> Hq 1 , the universe tends to be fully described 
by the de Sitter model. Here, we revisit the problem 
of quantizing the electromagnetic field in the de Sitter 
space as previously considered in Ref.[14j. We present 
an exact solution of the Schrodinger equation, construct 
the correspondent coherent states, and calculate various 
average values. The paper is organized as follows. In Sec. 
II we show how to obtain the correspondence between an 
electromagnetic field placed in a de Sitter background 
and a time-dependent harmonic oscillator. In Sec. Ill, 
we use the Lewis and Riesenfeld 12| method to obtain the 
solutions of the related Schrodinger equation, construct 
the coherent states, calculate the correlations between 
position and momentum and the uncertainty product. 
Section IV summarizes the results. 



DECOMPOSITION OF THE VECTOR FIELD 

Consider the electromagnetic field in a de Sitter space 
time. The action for the Gauge field is given by 



S = -\[d 4 x V^g^g^F^F^, 



(1) 



where the metric ds 2 = —dt 2 + e 2Hot dx 2 is diagonal and 
F^v = d^Av - d u A^, leading to 

g^g^F^F^ = - 25^e' 2H ° t F Ql F 0j + 
e -Wot S kl s iJ Fk . Fij _ 



Using the Gauge conditions 

A = 0, V-A = 
we obtain 



(2) 
(3) 



S= -- f Vxe^'C-e- 2 " *^) 2 + e-^&Aj) 2 ) 

(4) 
where the dot denotes the time derivative. Now we de- 
compose the field as 



A{x,t)= £ 



A=l,2 



d 3 k 

J2nf 



?(*) 



(5) 
where e represents a unit vector satisfying /c • e^ x ' = 
due to the Gauge condition and e* ' • e^ = 5 X A , and 
A = 1,2 stands for the two polarizations. By separating 
ujj, into real and imaginary parts 



(A) 
u k = - 



to 



X ' ^2) 



(6) 



the action reads 
5 = 



Ifdtf^Y, ~ eH0t [&* + fc2 ^] > ( 7 ) 



where i = 1 , 2 labels the real and imaginary parts of 



,W 



From Eq.fJT]) we obtain the Hamiltonian for the 



electromagnetic field 



-H t 



where 



H\ik — 



PXik 



'KPXik + & ^Aifcji 
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dq 



Xik 



- e Hot q Xlk , 



(8) 



(9) 



with p being the conjugate momentum. The classical 
equation of motion for the qth mode reads 



q\ik + H q\ lk + uj q Xlk = 



(10) 



Next, consider the classical harmonic oscillator with 
time-dependent mass and frequency given by the Hamil- 
tonian 



H(t) 



-\ + Jm(t) W 2 (t) 9 2 , 
m(t) 2 



(11) 



where [q,p] = ih. The equation of motion reads 

rn(i) . , . . 

q + ^q + u 2 (t g = o, 

m(t) 



(12) 



which is similar to Eq. (jlOjl if one considers that each 
mode of the electromagnetic field corresponds to the 
time-dependent harmonic oscillator with m(t) — 
and w = ke~ Hot . 



■,H t 



QUANTIZATION OF ELECTROMAGNETIC 
FIELDS IN THE DE SITTER SPACETIME 



Consider a time-dependent harmonic oscillator de- 
scribed by Eq. pip . It is well known that an invariant 
for Eq. (fTTj) is given by [l7| 



"3 



1 ) + (fiP - rap'g) 2 



(13) 



where q(t) sastisfy Eq. (IT2|) and p(i) satisfies the general- 
ized Milne-Pinney equation jl 51 ? , llq 



p + "f{t)f> + LJ 2 (t)p = 



1 



m 2 (t)p 3 
The invariant I(t) satisfies the equation 



(14) 



(15) 



and can be considered hermitian if we choose only the 
real solutions of Eq. ((H)) . Its eingenfunctions, (j> n (q,t), 
are assumed to form a complete orthonormal set with 
time-independent discret eigenvalues, A n . Thus 



I(j> n (q,t) = \ n (j) n {q,t) 



(16) 



with (</>„,</>„<) = S nn '. Taking he Schrodinger equation 

(SE) 



iK d -^=Hm(q,t) 



(17) 



where H(t) is given by Eq. (jTTJ) with p = —i%-S-. Lewis 
and Riesenfeld [12J showed that the solution ip n {q,t) of 
the SE (see Eq. ((T7|) ) is related to the functions <f> n (q,t) 
by 

lM?>*) = e*" {t) tf»(?>*) (18) 

where the phase functions n (t) satisfy the equation 



Mi) , 

n -r = \4>n{q,t) 



tol-H® 



4> n {q,t)). (19) 



The general solution of the SE may be written as 



i>n(q,t) = J2 C n ei6n(t)( t>n(<l,t) 



(20) 



where [a^(t),a(t)] = 1. In terms of a(t) and aJ(t) the 
invariant I (see Eq. (jT3l)) can be written as 



I = h{a{t)a\t) + 1/2). 



(26) 



Let \n, t) be the eigenstates of /. Therefore the follow- 
ing relations hold 



where c„ are time-independent coefficients. Now, using 
an unitary transformation and following the steps as in 
Ref. [13 we find 



Mq,t) =e i9n(t) 



where 



7r 1 /2n 1 / 2 n!2> 



1/2 



exp 



im(t) 
2ft 



2ft m(t)p 2 (t) 



q 2 x (21) 



" '.*? 



a(t)\n, t) =y/n\n — 1,4), 



a + (i)|n,i) =V(n+l)|n+l,t), 
I|n,t) =ft(n + l/2)|n,t). 



(27) 
(28) 
(29) 



Since the coherent states for / can be easily constructed, 
the coherent states for the time-dependent harmonic os- 
cillator are straightforwardly obtained: 

2 OO 

MHe-^^T-^e^M), (30) 

n— 

where n is given by (|22p . and the complex number a(t) 
satisfies the eigenvalue equation 



»(*)="(»+ o) 



2'J to m(t')p 



dt' 



(22) 



A n = (n+i)ft, and -ff„is the Hermite polynomial of order 
n. 

For m(t) = e Hot and w(4) = ke~ Hat the auxiliary 
Equation (TT4")) reads 



-2H t 



p + H oP + k z e-' nat p 



(23) 



Following the procedure described in Ref. [151 ? , [iff , we 
find p = 1/fc 5 , and from Eqs. (|22|) and (|2"Tj) we have 



ipn(q,t) 



k 



exp 



7rftn! 2 2™ 
kq 2 



* e -T%("+!)(e- H " t0 -e- Hot ). 



2/i 



i?i, 



ft 9 , 



(24) 



which, except for the phase factor, is similar to the well- 
known wave function for the time-independent harmonic 
oscillator. 

The coherent states for the time-dependent harmonic 
oscillator (see Eq. ([IT])) are constructed as follows [18J . 
Consider the time-dependent creation a'(t) and annihi- 
lation a(t) operators defined as 



ot(t) 

a(t) 



1 



-i(pp-m{t)pq) 

+ i(pp - m(t)pq) 



(25) 



with 



and 



a(t)\a, t) — a(t)\a,t), 



a(t) = a{t )e 2l9o{t) 



0o(t) 



dt' 



2 J t B m(t')p 2 (t')- 



(31) 



(32) 



(33) 



The fluctuations in q (Aq), p (Ap), and in the uncer- 
tainty product AqAp of the coherent state \a,t) are 



(34) 



A 1 a = \I^P, 



ftl 



and 



A P , = \/~(l + m 2 (i)pV) 1/2 , 



Aq a Ap a = - (1 + m 2 (t)p 2 p 2 ) 1 / 2 , 



(35) 



(36) 



respectively. 

For p — 1/k 1 / 2 , we have m{t)pp — and as a conse- 
quence Aq a Ap a is minimum (ft/2). In this case, the 
coherent states \a,t) are "true" minimum-uncertainty 
states. Otherwise, if m(t)pp ^ 0, the states \a,t) are 



not minimum-uncertainty states and are equivalent to 
the well-known squeezed states, as pointed out in Ref. 
19]. This is an interesting result since the minimum 
uncertainty product is assumed to be satisfied only for 
time-independent harmonic oscillator, unless the solution 
of Eq. (fT4l is a constant. 

Next, let us analyze the (q) a and (p) Q time evolution 
and its phase diagram. By setting a(to) = u + iv we find 



A-7o 



2hp [u cos(20(i)) - v sin(26» (t))] (37) 



and 



Ap a 



— h ump cos(28o(t)) + 
P 



vmp ) sin(20o(t)) 



(38) 



The constants u and v are determined from the initial 
conditions (q(t )) a = q and (p(t )) a = p = m(t )v . 
Then we find 



and 



(Q) a 


—qocos 


(e~ Hot -e~ H ° to ) 

Ho 






Po ■ 
-—sin 
k 


' k 
ft ( e -H t e -H to\ 




(P) a 


=pocos 


k 

_( e -H a t _ e -H to\ 

Ho 






kqosin 


— (e~ Hot -e' Hoto ) 
Ho K J 





(39) 



(40) 



Figures HJa) and [TJb) show the time dependence of 
(q) a and (p) a . Since we want to work in a de Sitter uni- 
verse, we used the values: to ~ -f^o" 1 — 10 18 , k = 10~ 8 , 
go = 1-0, and po = 1.0. From Fig. QJa) one can see that 
the system oscillates back and forth between the classi- 
cal turning points with an increasing period and constant 
amplitude. For t sa 2AH , q and p reach their asymptotic 
values given in Eqs. (|4*T|) and (|4*2"j).The phase diagram 
is shown in Fig. [He). Even though the dissapative na- 
ture of the system (the total energy (E ? 



-H t 



_ ot )),_it 

behaves like the usual time-independent harmonic oscil- 
lator (E =constant). This can be seen from the relation 
B = y/2E/X, where B is the amplitude of motion and A 
the spring constant. Since A = k 2 e~ Hot and E ~ e~ Hot , 
B is a constant. As t increases, the oscillation frequency 
u(t) decreases while the oscillation period increases lead- 
ing to the behavior shown in Fig. 1(a). It is also worth 
noticing that for t >> (1/Hq), the oscillations disap- 
pear and the motion of both (q) a and (p) a tends to be 
a constant. This kind of behavior is expected in an eter- 
nally expanding universe. Therefore, in the above figure 
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FIG. 1. (a) The time evolution of (q), (b) Time evolution of 
(p), and (c) The phase space. We considered k = 1CP 8 in all 
panels 



we consider our time in the range to = 10 18 < t to il- 
lustrate the time evolution in this model. As said in the 
introduction, for the ACDM model, our universe gets 
closer and closer to de Sitter, this justifies the inclusion 
of large times (t —t oo) in the figure. In Fig. [2] we show 
the behavior of the generalized coordinates and the phase 
space for a smaller value of the wave vector k = 10~ 4 . 
From this we can see that the modes with shorter wave- 
length oscillates for a longer time. This is expected since 
as the universe expand, the wavelengths are stretched 
and eventually disappear. 

In the limit when t — > oo, (q) and (p) are expressed 
as 



(?) a = c os ( — < 



-H()t( S 



-Hoto 



Hn 



(41) 



and 



(p) a = cos [ —e 



— ksin I — —e 
Ho 



-H®to 



(42) 



From the expressions we see that the coordinates depends 
explicitly on k. We show in Fig[3]the behavior of the sys- 
tem when t goes to infinity. The (q) and (p) dependence 
on k is shown in panels a and b of the figure, while the 
phase space is shown in Fig. [3jc). In this Figure and the 
above expressions we see that for large values of k the (p) 
axis is more stretched while for small values this is so for 
the (q) axis. 

Pedrosa et al. [19( have combined linear invariants 
and the LR method to obtain the exact wave function 
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FIG. 2. (a) Time evolution of (g)(b) Time evolution of (p)and 
(c) Phase space for k — 1CP 4 
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FIG. 3. For t — s- oo : (a)fc dependence of (q) (b)fc dependence 
of (p) and (c) k dependence of phase space. 



for a particle trapped by oscillating fields, written in 
terms of Mathieu functions. They calculated AqAp and 
the quantum correlation between q and p, defined by 
c i,i = \{{lP + Pi)) ~ (<z)(p)- They are related through 
the equation 



AqAp 




(43) 



which shows that AqAp is minimum whenever ci,i = 0, 
as it happens for c\ t \ calculated in the coherent state 
\a,t) 1 i.e., (ci j i) a = 0. The fact that ci.i = does not 
mean that q and p are uncorrelated. In order to verify the 
correlation between q and p one should study the function 
c n , m = \{{q n p m +p m q n )) - (q n )(p m )- For the coherent 
state | a, t), we find that (02,2)0 = — h /2 , indicating 
that q and p even assumed as "classical" quantities are 
correlated. 

The uncertainty product and correlations in the 
state ip n are more easily calculated using the relation 
|V'n( < 7^)) = e l9n ^\n,t). They are given by 



Aq^ n Ap^ 



(ci,i)v-„ = 0, 



and 



( c 2,2), 



n + n + — ) h 



(44) 



(45) 



(46) 



case ?/>o is given by 



'■"" / - /, = ^) 4e " aS ° e " H ° lexp (-^) ffo \fn q , 



(47) 



in the coordinate representation of the coherent state. 



CONCLUDING REMARKS 

In this paper we used the Lewis and Riesenfeld method 
to obtain the time-dependent Schrodinger states emerg- 
ing from the quantization of the electromagnetic field in 
the de Sitter spacetime. It is well-known that a challenge 
in obtaining the exact solution (see Eq. ([21~]) ) for the SE 
with H given in Eq. (|L 1 [) . is the solution of the Milne- 
Pinney (in terms of the universal scale factor of universe 
a(t)) equation 



-p + uj 2 (t)p 



a 2 (t)p 3 



We noticed that Eq. (|37p is satisfied for n — 0. In this 



For a given a(t), different eras of universe can be in- 
vestigated. For the de Sitter universe a(i) = e Hot , and 
P=^-pii and the solution for ty n corresponds, except for 
the phase factor, to the well known wave functions of the 
time-independent harmonic oscillator. 

We have constructed the "true" coherent states, \a,t), 
whose coordinate representation is given by Eq. (l47l) . We 
verified that Eq. (|4"3")) holds for \a,t) . We calculated the 
quantum fluctuations in the coordinate and momentum 
as well as the quantum correlations between the coordi- 
nate and momentum in the state ^ n (q, t). We calculated 



the quantum fluctuations in the coordinate and momen- 
tum as well as the quantum correlations between the co- 
ordinate and momentum in the state ^ n (q,t). 

We analyzed the time behavior of q and p, as well as 
the phase diagram q and p (see FigsOJa)-(c) and[5fa)-(c) 
). The phase diagram indicates that the properties of the 
light propagation through the de Sitter spacetime can be 
analyzed in terms of the well known time-independent 
harmonic oscillator in real time. 

At last, we would like to point out that the procedure 
described can be used to trace the present properties of 
the quantum electromagnetic field back to the recom- 
bination era. This would be a much more interesting 
phenomenological result. We are currently investigating 
the properties of the quantum electromagnetic field for 
a(t) ~ fi and a(t) ~ t~s to give a more complete de- 
scription of the cosmological evolution of the quantized 
electromagnetic field. 
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